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Abstract. There is a high demand of space-efficient algorithms in built- 
in or embedded softwares. In this paper, wc consider the problem of 
designing space-efficient algorithms for computing the maximum area 
empty rectangle (MER) among a set of points inside a rectangular re- 
gion TZ in 2D. We first propose an inplace algorithm for computing the 
priority search tree with a set of n points in TZ using O(logn) extra bit 
space in 0(ri log n) time. It supports all the standard queries on prior- 
ity search tree in 0(log^ n) time. We also show an application of this 
algorithm in computing the largest empty axis-parallel rectangle. Our 
proposed algorithm needs 0{nlog^ n + m) time and O(logn) work-space 
apart from the array used for storing n input points. Here rn is the 
number of maximal empty rectangles present in TZ. Finally, wc con.sider 
the problem of locating the ma^dmum area empty rectangle of arbitrary 
orientation among a set of n points, and propose an O(n^logn) time 
in-place algorithm for that problem. 



1 Introduction 



Though memory is getting cheap day by day, still there are massive demand 
for the low memory algorithms for practical problems which need to be run on 
tiny devices, for example, sensors, GPS systems, mobile hand-sets, small robots, 
etc. Also, now-a-days the data available in several problems itself is huge. So, 
the practical algorithms for the data-streaming or data-mining problems must 
be space-efficient. For all these reasons, designing low-memory algorithms for 
practical problems have now become a challenging task in the algorithm research. 

In computational geometry, designing the in-place algorithms are studied only for 
a very few problems. For convex hull problem in both 2D and 3D, space efficient 
algorithms are available. In 2D, the best known result is a 0{n\ogh) algorithm 
with 0(1) extra space [5]. Bronnimann et al. [4] also showed that the upper 
hull of a set of n points in 3D can be computed in 0(n log^ n) time using 0(1) 
extra space. In the same paper it is also shown that for a parameter s satisfying 
clog^ n < s < n, (c > is a constant), if O(^) space is allowed, then the convex 
hull for a set of n points in 3D can be computed in in 0{ns) time. Vahrenhold 



[14] proposed an 0{n^ logn) time and 0(1) extra space algorithm for the Klee's 
measure problem, where the objective is to compute the union of n axis-parallel 
rectangles of arbitrary sizes. Bose ct al. [3] used in-place divide and conquer 
technique to solve the following three problems in 2D using 0(1) extra space: 
(i) a deterministic algorithm for the closest pair problem in O(nlogn) time, (ii) 
a randomized algorithm for the bichromatic closest pair problem in O(nlogn) 
expected time, and (iii) a deterministic 0(rilogn + k) time algorithm for the 
orthogonal line segment intersection computation problem. For arbitrary line 
segments intersection computation problem, two algorithms are available in [8]. 
If the input array can be used for storing intermediate results, then the problem 
can be solved in 0((n + fc)logn) time and 0(1) space, but, if the input array 
is not allowed to be destroyed, then the time complexity increases by a factor 
of logn, and it also requires O(log^n) extra space. Regarding the empty space 
recognition problem, it needs to be mentioned that all the Delauney triangles 
among a set of n points can be computed in O(n^) time using 0(1) space [2]. 
This, in turn, recognizes the largest empty circle among a point set with the 
same time complexity. 

In this paper, we propose an in-place algorithm for constructing a priority search 
tree T [10] with the set of points P in 7?., jPj = n. This needs an additional 
O(logn) bits. We show that, the standard queries on priority search tree can 
be performed in it in O(log^n) time. Priority search tree is a very important 
paradigm in geometric algorithms, as it has several important applications. Thus, 
our algorithm may aid in several problems in memory- restricted environment. 

As an immediate application of our inplace priority search tree, we have consid- 
ered the computation of largest empty axis-parallel rectangles among the points 
in P. Our proposed algorithm runs in 0(TO-|-nlog^ n) time using O(logn) extra 
space. Here m is the number of maximal empty axis-parallel rectangles (MERs) 
in TZ. By maximal empty axis-parallel rectangle (MER), we mean an empty 
axis-parallel rectangle that is not containined in any other empty axis-parallel 
rectangle. 

The problem of computing the largest MER was first introduced by Namaad et 
al. [11]. They showed that the number of MERs' (m) among a set of n points 
may be O(n^) in the worst case; but if the points are randomly placed, then the 
expected value of m is O(nlogn). In the same paper, an algorithm for identifying 
the largest MER was proposed. The worst case time complexity of that algorithm 
is 0(TOm(n^,TOlogn)). Orlowski [13] proposed an easy to implement algorithm 
for finding the largest MER that runs in 0(m + nlogn) time. It inspects all 
the MERs' present on the floor, and identifies the largest one. The best known 
algorithm for this problem runs in 0(n log^ n) time in the worst case [1]. Same 
time complexity result holds for the recognition of the largest MER among a set 
of arbitrary polygonal obstacles [12]. Recently, Boland and Urrutia [6] gave an 
0(n logn) time algorithm for finding the largest MER inside an n-sided simple 
polygon. All these algorithms use 0(n) extra space. 



Finally, we considered the problem of designing an in-place algorithm for com- 
puting the largest empty rectangle of arbitrary orientation among a set of n 
points. It takes O(n^logn) time with an 0(1) additional workspace. The best 
known algorithm for this problem in the literature runs in 0{n^) time and O(n^) 
space [7]. 

2 In-place priority seeirch tree 

Let P = {pi,P2, ■ ■ ■ ,Pn} be the given set of points in 2D, where 2*^"^ < n < 
2''. The array that stores the points in P, is also referred to as P. The i-th 

array location will be referred as P[i]. We now define the priority search tree T 
recursively in a way that suits our in-place implementation. 

The tree T has exactly k levels. The level of root in T is considered to be the fc-th 
level. The root represents the entire set of points P, and it stores the point p* , 
where y{p*) = maxpg p y{p). Its two children are the priority search tree with the 
set of points and P^, where P^ and Pr are defined as follows: (i) let m = 2*^"^. 
Compute m-th order statistics X(^m) among the a;-coordinates of the points of P. 
The set P^ = {p e P I p 7^ P*,Xp < x^„i)}, and Pr = {p E P \ p ^ p* ,Xp > 
Thus, \Pe\ = 2*^-1 - 1, |P^| = n - 2'=-\ and P^ U P^ = P \ {p*}. In our modified 
definition of priority search tree, we assume that each node at level i {i ^ 0) 
represents a tree of size 2* — 1 except the rightmost node in that level. Observe 
that, at the level k — 1, each of the two subtrees of Pi are of size 2*^"^ — 1. But 
for Pr, we may not always be able to construct two subtrees if \Pr\ < 2*^~^. 
In that case, it has only the left subtree rooted at the point having maximum 
2/-coordinate in the point set Pr', otherwise it has two subtrees. 

In general, at any node of the z-th level, the root is as defined earlier. If P' 
denotes the set of points represented by that node, then the number of children 
of that node is (i) zero or (ii) one or (iii) two depending on whether (i) |P'| = 1, 
or (ii) 1 < |P'| < 2*-i or (in) |P'| > 2'"^. In Case (i), it has no subtree. In Case 
(ii), it has only the left subtree, and its size is [P'[ — 1. In Case (iii), it has both 
left and right subtrees, and their sizes are 2*~^ — 1 and |P'| — 2*~^ respectively. 
We maintain an array TAG with 2 log n cells, indexed by the levels of the tree T. 
Each cell is of size 2 bits. The TAG[i] is set to or 1 or 2 depending on whether 
Oi = the number of nodes at level i of T is equal to 26'i_i or 26'j_i — 1 or 
20i-i — 2. Note that, unlike the usual priority search tree [10], at each node the 
discriminating x-value among the points in two subtrees is not stored. Here, the 
method of deciding the search direction from a node will be decided by observing 
the inorder predecessor and successor of that node in T. 

2.1 In-place organization of T 

T is organized in a heap like structure. All its leaves appear in the same level; 
but unlike heap, a non-leaf node of T may have zero or one or two child(ren). 



In a particular level i, at most one node may have less than two children, and 
if such a node exists, it is the right-most node in that level. The root of T is 
stored in P[l]; it has always two children, stored in P[2] and P[3] respectively. 
In general, if all the nodes in T have two children except the leaves, then the 
children of P[j] reside at P[2j] and P[2j + 1]. But, since at most one node in 
a level of T is permitted to have less than two children, we use TAG bits to 
compute the address of the children of a node. If a node at level i, and stored at 
P\j], has two children, they are available at P[2j — Y^'^^^_^_l TAG[a]] and P[2j — 

Yla=i+i TAG[a] + 1] respectively. If P[j] corresponds to the right-most node at 

level i, and TAG[i\ = 1, then its only child resides at P[2j - J2t=i+i TAG[a\]. 

2.2 Creation of T 

We first initialize r^G[i] = for all i = l,2,...,k, where k = [lognj. The 
computation of 7~ is done in a breadth-first manner. In other words, all the nodes 
in a particular level i are computed prior to computing the nodes of level for 
alH = 1, 2, . . . , fc. We compute p* = P[j] (say), where y{P[j]) = maxf^i y{P[i]) 
as the root of T, and store it in by swapping P[l] and P\j]. Next wo sort 
P \ {p*} in an in-place manner with 0{n) data movement [9]. Its first 2*'^^ — 1 
elements form the set Pi for the left subtree Te, and the remaining n — 2^^^ — 1 
points form the set Pr for the right subtree %■■ Next, we identify p*^ and p*, the 
roots of Ti and %. (as defined for T). The point p*^ (resp. p*) is swapped with P]2\ 
(rcsp. -P[3]). Again we sort the points in P \ to compute the nodes 

in the next level. Note that, at this level, both the children of p| exists; but the 
number of children of p* may be zero or one or two. If the number of children 
of p* is one or zero, r^G[2] is set to 1 or 2. The same process continues up to 
the fc-th level. Since, at each level, a sorting is involved, we have the following 
result: 

Lemma 1. For a given set P of n points, the tree T can he constructed in 
O(nlog^n) time. 

2.3 Traversal of T 

The traversal in T starts from its root (at P[l]). At each step, it moves towards 
one of its children. If T is full, the children of a node (point) stored at P[j] 
(at level i of T) are available at P[2j] and P[2j -|- 1]. But, since T may not be 
full, we need to use TAG bits attached to each level of T for the traversal. We 
maintain an integer location A during the traversal of T- While moving from 
level i to level i + I, we add TAG[i] with A. It is already mentioned in the 
earlier subsection that, the left (resp. right) child of the node P[j] are available 
at location 2j — A (rcsp. 2j + 1 — A). Again, if P[j] is the right-most node of a 
level of T, it may have zero, one or two children, and this information is available 
at TAG[i]. If P[j] is the right-most in its level i, and it has only one child, then 



the algorithm has to move towards its left child irrespective of the requirement 
(of moving towards left or right) in the algorithm. 

3 Standeird queries on priority seeirch tree 

We now show that the standard queries on the priority search tree [10] can be 
performed in T also in 0(log^ n) time with 0(1) additional space. 

3.1 MinXInRectangle(a;o, cci, j/i) 

Three real numbers xq, xi and jji arc given. The objective is to find the point 
p* = {x* ,y*) e P with minimum x-coordinate among those points p = {x,y) € P 
satisfying xo < x < xi and y > yi- 

Since the a;-coordinate of the partitioning line at each node of T is not stored 
as in [10], the search direction from a node p S T is decided by its inorder 
predecessor p~ and inorder successor p+ . 

While executing this query with the interval [a;o,a:i], we first find the discrimi- 
nant node IT € Tsuch that x{-K~) < xi and x{n^) > xq. If y(7r) < yi, then report 
that the search can not output a feasible point satisfying the query; otherwise 
the search proceeds to answer the query. We initialize two locations p* and A* 
with TT and A, where p* will contain the final answer, and A is the sum of TAG 
bits computed during the traversal up to the node tt. The search proceeds in 
both the subtrees of tt. The traversal in the left subtree of n starts with p = w. 
The actions taken at each node p on the search path, and the choice of its child 
for the next move is decided as follows. 

• If y{p) < yi, the search stops. Otherwise, execute the following three steps. 

• If Xq < x(p) < Xi then assign p* ~ p 

• If the x{p~) < Xo, then set p = right child of p. 

• If the x{p'^) > Xo, then set p = left child of p 

Using a similar procedure we traverse the right subtree of tt starting with p = i:, 
and restoring the value of A at node tt, which is stored in A*. Finally, p* is 
reported as the answer to the query. 

Time complexity: Since here the search direction from a node q is guided by 
x{q~) and x{q~^), each move from a node to its descendant in the direction of 
the search takes O(logn) time. Note that, while, searching q~ or g+ of a node 
q, we copy Z\ of g at a scalar location A' , and perform the search as mentioned 
in subsection 2.3. 

Since the total number of nodes to be traversed to report the answer or the 
non-existence of a feasible solution is O(logn), we have the following lemma: 



Lemma 2. Using the in-place priority search tree on a set of n points, the 
]V[inXInRectangle(a;o,a;i,yi^ query can he answered in O(log^n) time using 
0(1) extra space. 



3.2 MaxXInRectangle(a;o, a;i, yi) 

Three real numbers xq, x\ and yi arc given. The objective is to find the point 
p* = {x*,y*) € P with maximum x-coordinate among those points p = (x.y) e 
P satisfying Xq < x < Xi and y > yi- This query can be answered in a similar 
manner as in MinXInRectangle query with same time and space complexity. 



3.3 MaxYInXRange(a;o,a;i) 

Given a pair of real numbers xq and xi, find a point p* = {x*,y*) whose y 
coordinate is maximum among all the points in P satisfying xq < x < Xi. Here 
the algorithm is essentially the same as in MinXInRectangle query. Here if a 
node satisfies xq < x < xi during the search for the discriminant node, the search 
stops reporting that point. Otherwise, we need to search separately both the 
subtrees of the discriminant node till a point is obtained satisfying xq < x < xi. 
Surely, time and space complexities of MinXInRectangle query hold here also. 

3.4 EnumerateRectangle(a;o, a^i, yi) 

Three real numbers xq, x\ and y\ arc given. The objective is to identify all the 
points p = {x,y) £ P satisfying xq < x < xi and y > yi- Here, the discriminant 
point TT is found as in MinXInRectangle query. In this path, if there exists any 
point satisfying the given constraint, then report it. Next, perform an inorder 
traversal in the subtree rooted at tt to identify all the points satisfying the desired 
condition. During the inorder traversal, (i) if a node with inorder predecessor 
having x-coordinate less than xq is reached, its left subtree is not traversed, 
similarly (ii) if a node with inorder successor having x-coordinatc greater than 
Xi is reached, its right subtree is not traversed, and (iii) if a node is reached 
whose 2/-coordinate is less than yi, then its both the subtrees are not traversed. 

Also note that, from a node P[j] at level i, the index of its parent (at level 
in the array P is computed as J _ TAG[i — 1]. At each movement 

from a node at level i to its parent at level i + 1 in T, we need to update A by 
A — TAG[i]. Thus we have the complexity results of this query as stated below: 

Lemma 3. Using the in-place priority search tree on a set of n points, the 
EnumerateRectangle(lro, a;i, j/ij query can be answered in 0{m + \og^ n) time 
using 0(1) extra space, where m is the size of the reported answer. 



4 Lcirgest empty rectangle 



We now concentrate on our main problem of computing the largest empty axis- 
parallel rectangle among the point set P stored in an axis-parallel rectangular 
region TZ. The algorithm consists of two phases: top-down and bottom-up. Each 
phase consists of n passes. In each pass of the top-down phase, we identify all the 
MERs with the point stored at the root node of T, on its top boundary, and then 
delete the root from T. Thus, a new point having maximum y-coordinate among 
the remaining points in P becomes the root. The same algorithm is repeated to 
compute MERs with the new root at their top boundary. The deletion of the 
root of T is explained in detail. After the deletion, a space in T becomes empty. 
Actually, we store the deleted root of T in that location. We show that, it does 
not affect the correctness of our algorithm. Thus, after the completion of the 
top-down phase, all the points in P are present in the array P, and we can 
execute the bottom-up phase with all the points in P stored in the same array. 
The bottom-up phase is exactly similar to the top-down phase. Here, in each 
pass, all the MERs with bottom boundary passing through the point stored in 
the root node of T are identified. We now explain the top-down phase in detail. 

4.1 Top-down phase 

We now explain the processing of the root p* = P[l] of T. Let Xmin and x^ax be 
the left and right boundary of TZ respectively. We use two double-ended queues 
Qi and Qr to store the points encountered in the two sides of p* during the 
traversal of T. It stores some already visited points of T for the future processing, 
and will be clear from subsequent discussions. The insertion and deletion in both 
the queues can be performed in both of their ends. At the begining of each pass, 
Qi and Qr are empty. We define a curtain with horizontal span I = [xmin, Xmax]', 
its top boundary is fixed at p*. Let pi and Pr be the two children of p*. If pe and 
Pr are in different sides of p* , then both the points are pushed in their respective 
queues. But, if and pr are in the same sides (say left) of p*, then two situations 
need to be considered: (i) if y{pe) > y{Pr), then both pr and pe are pushed in 
Qi in this order. Otherwise, pr is pushed in Q^, and pf is ignored. 

Next time onwards, the point p for the processing is the one having higher y- 
coordinate among the front elements of Qe and Qr- While choosing the point p 

for processing, it is first deleted from the respective queue, and then the MER 
is reported as stated below. It also causes updating of the queues Qi and Qr- 

4.2 Processing the topmost queue element 

Let X = [a, f3] denote the horizontal span of the curtain. We first report an MER 
with the horizontal span I, and vertical span [y{p*),y{p)]- I is truncated at p, 
so that p* lies in the updated I. Here also, we will use pe and Pr to denote the 



Fig. 1. Three difFerent cases while processing a point p in the top-down phase 



children of p. Depending on the position of pe and Pr with respect to p* and p, 
one or both of pi and Pr are put in the appropriate queue as described below. 

Without loss of generality, let us assume that p is to the left of p* , and pi and 
Pr be its two children. Here, three cases may arise: (i) both p^ and Pr are to the 
left of p (Figure 1(a)), (ii) both pi and Pr are to the right of p (Figure 1(b)), and 
(iii) Pi and Pr are in different sides of p (Figure 1(c)). The actions taken in the 
three different cases are stated below. 



- paitition line for p and q 
partition line for pf ami p,. 



(a) (b) 

Fig. 2. Demonstration of Case (i) 



Case (i) Here, both pi and Pr are not inserted in Qi. However, we need to 
follow the right links starting from pr until (a) a point p' is found inside X, 
or (b) the bottom boundary of the floor is reached (i.e. the index of the right 
child of the current node on the traversal path is greater than the size n of 
the array P). In case (a), surely we have x{p') < x{q) where q is the element 
at the front-end of Qg^. Now, if y{p') > y{q), we insert p' at the front-end 
of Qi (see Figure 2(a)). But, if y{p') < y{q), we ignore p' , or in other words, 
do not insert it in Qg (see Figure 2(b)). 

Case (ii) Here three different situations may arise: (a) both pi and Pr are to 
the left of p* , (b) both pi and Pr are to the right of p* , and (c) pi and Pr 
are in different sides of p*. In subcase (a), if y(pi) > y{Pr), we insert pr 
and pi at the front-end of Qg in this order (see Figure 3(a)); otherwise only 
Pr is put in Qi, and pi is ignored. In subcase (b), we insert both pr and 
pi at the rear-end of Qr- This may need deleting elements of Qr from its 

^ The reason is that tlie point is p' is in the same partition of p and the point q is 
entered in Qi by the right sibling of p or the right child of some ancestor of p. 



Fig. 3. Demonstration of Case (ii) 



rear-end. Such a situation is demonstrated in Figure 3(b). Here a and b are 
already present in Qr] while inserting pr, b is deleted from its rear-end since 
y{b) < y{pr)- Next, is also inserted in Qr in the same manner. Note that, 
while inserting p^, pr may also be deleted (see Figure 3(c)). In subcase (c), 
Pi and Pr will need to be inserted in and Qr respectively. Note that, if 
is non-empty, then pi appears to the left of all the elements present in 
Qi. Thus, if y{pi) is greater than the y-coordinate of the element present at 
the front-end of Q^, then pi is inserted at the front-end of Qi. Similarly, if 
there is any element in Qr^ Pr will be left to all of them. Here, Pr must be 
added at the rear-end of Qr ; if needed, some element of Qr are deleted from 
its rear-end. The tiny points in Figure 3(d) are already present in the Qr- 
Note that, if this situation arises, at least one of Q^ and Qr must be empty. 
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Fig. 4. Demonstration of Case (iii) 



Case (iii) Here, first pr is inserted in either at the front-end of Qe or at the 
rear-end of Qr depending on whether pr is to the left or right of p* (see 
Figure 4). Note that, here pg must be ignored. But as in case (i), we need 
to proceed following the right links starting from pi to reach either a point 
p' inside I or the bottom boundary of the floor. If y{p') > y{q), where q is 
at the front-end of Qg, then p' is entered in Qg at its front-end (see Figure 
4(a)); otherwise p' is ignored (see Figure 4(b)). Figure 4(c), shows a situation 
where Pr is to be inserted at the rear-end of Qr- 

If p has a single child, it is inserted in the queue in a similar manner. If p has no 
child, no special action needs to be taken; we only proceed to process the next 
point in the queue. 



The similar set of actions are taken for processing the point p when it is at the 
right side of p* . The execution continues until both the queues become empty. 
Then the last (only one) MER is reported with the resulting curtain I as the 
horizontal span, and vertical span defined by p* and the bottom boundary of TZ. 

4.3 Deletion of the root of T 

After the completion of a pass, we reorganize the tree T by deleting the point 
at its root as follows, without computing it afresh. 

We start from the root by assigning i = 1. At each move we consider both the 
children of P[i], and choose the one, say P[j] having maximum y-coordinate. 
The tie, if arises, is broken arbitrarily. If P[i] has only one child, its index is 
taken in j. We swap P[i] and P[j]. and the algorithm proceeds by copying the 
value of j in i. Finally, when no child of P[i] is found, the algorithm terminates. 
Here the following facts need to be mentioned: 

Fact 1 The point at the root is logically deleted, but it still remains in the array 
P. This is essential, since we need to execute a bottom-up phase with all the 
points in P after the execution of the top-down phase. 

Fact 2 Usually, while traversing along a path ofT, the bottom of the floor TZ is 
recognized, when a leaf is reached, or in other words, the index of the child of the 
said node along the desired direction is greater than the array size n. But, such 
a method may fail from the second pass onwards due to our scheme of deleting 
the root. Here apart from the usual way, the leaf is also detected if a node (point) 
is reached whose y-coordinate is greater than the y-coordinate of the root ofT. 

Fact 3 During the deletion of the root ofT, when a point is moved from a child 

node to its parent, it still remains in its own partition according to the partition 
line defined at that node at the time of creation of the T. Thus, binary search 
property according to the x-coordinate values in the present T still remains valid. 
Moreover, the point having m,aximum y-coordinate in a partition is at the root 
of its corresponding subtree. Thus the modified T is still a priority search tree. 

Fact 4 The tree T may no longer remain balanced after the deletion of some 
points from T. Or in other words, leaf node may appear in different level. How- 
ever, the length of each search path will still be bounded by O(logn). 

By the virtue of Facts 1, 2, 3, we can execute the subsequent passes to identify 
all the MERs with top boundary aligned at the point staying at the root node 
of T in that pass. Fact 4 says that the time required for deletion of root at the 
end of each pass is O(logn). 

4.4 Complexity Analysis 

Lemma 4. A single pass of the top down phase needs 0(/LI + logn) time in the 
worst case, where is the number of reported answers in that pass. 



Proof. In each pass, when an MER is reported, at most two points are inserted 
in the queue. Thus, the number of points inserted in the queue is at most 2/i. 
Though the insertion of a point in the front-end always takes 0(1) time, the 
insertion of a point in the rear-end may need some deletions prior to that. But, 
since the total number of deletions in a pass is at most equal to the number of 
insertions in that pass, and no point is inserted twice in a pass, the time required 
for a single pass in the top-down phase is 0(/x). The time for the deletion of the 
root at the end of a pass needs O(logn) time. Thus, the result follows. □ 

Lemma 5. \Qi\ and \Qr\ can be at most O(logn) at any instant of time. 

Proof. The result follows from the fact that at any instant of time, Qe (resp. 
Qr) contains at most two points of a particular level in T. We justify this claim 
assuming the contradiction. Let 11 {\n\ > 2) be the set of points at the same 
level of T that are present in Qe at an instant of time. Surely, the parents for 
all of them are not the same, and they are inserted when their parents produced 
MERs. If we consider the x-coordinates of these points as well as their parents, 
there may exist at most two points (more specifically, the right-most two points) 
in 77, say tti and 7r2, such that x^tti) > .x(7r) and x{Tr2) > x^tt), where tt is the 
right-most one among the already processed parent nodes. Moreover, if such a 
situation arises, then tti and -1T2 are the children of tt. Thus, the other points of 
n have a;-coordinate less than x{7r), and hence they can not belong in Qg. □ 

Theorem 1. The time complexity of our algorithm for identifying the maximum 
area/perimeter axis-parallel rectangle among a set ofn points is 0{rn+n\o^ n), 
and it uses O(logn) work-space apart from the array P containing input points. 

Note: The time complexity of the algorithm can be reduced to 0(m + nlogn) 
if T can be constructed for the first time in O (nlogn) time avoiding the sorting 
at every level. 

5 Finding MER of cirbitrary orientation 

We now propose an in-place algorithm for finding maximum area empty rect- 
angle of arbitrary orientation among a set of points in P. The problem was 
addressed by Chaudhuri et al. [7] . They introduced the concept of PMER; it is 
the maximum area empty rectangle of any arbitrary orientation whose four sides 
are bounded by the members of P. It is shown that the number of PMERs is 
bounded by 0{n^) in the worst case. It follows from the following observation: 

Observation 1 [7] At least one side of a PMER must contain two points from 
P, and other three sides either contain at least one point of P or the boundary of 
TZ. A corner incident at the boundary of R implies that the corresponding sides 

contain that boundary point. . 

The worst case time complexity of the algorithm proposed in [7] is O(n^), and 
it takes O(n^) work-space for executing the algorithm. Our algorithm uses 0(1) 
work-space but its worst case time complexity is 0(n^ log n). 



5.1 Algorithm 



Observation 1 plays the central role in our algorithm. We consider each pair 
of points Pi,Pj S P, and compute all the PMERs with one boundary passing 
through Pi and pj. We assume that the points in P are increasingly ordered with 
respect to their x-coordinates; if tie occurs, then those points are increasingly 
ordered with respect to their y-coordinates. Two variables i and j are used to 
indicate the pair of points chosen for the processing. We choose different values 
of (i, j), i = 1, 2, . . . , n — 1 and j = i + 1, i + 2, . . . , n in this order. Each time the 
pair (i,j) is chosen, pi and pj are swapped with pi and P2 respectively. 

We execute the procedure Process{pi, pj) to compute all the PMERs with 
their one boundary passing through {pi,pj). Note that, after the execution of 
Process the points in P will not be in the increasing order of their co- 
ordinates as mentioned above. Thus, in order to choose the next pair for the 
processing, we need to sort the array P again with respect to their coordinates. 

Process (pi, pj): Consider the straight line £ij passing through Pi,Pj. It is trun- 
cated by the boundary of TZ at its two ends. The points Pi and pj are assumed 

to be stored in P[l] and P[2] respectively; the other points are split into two 
subsets according to the side of £ij it belongs. If Pi and P2 be the sets of points 
that are below and above lij respectively (|Pi| -|- IP2I = n — 2), then the points 
in Pi are stored in P[j],j = 3, . . . , |Pi| + 2, and the points of P2 are stored in 
P[j],j = |Pi| +3, . . . ,n. We use the following procedure to partition P\{pi,pj} 
into Pi and P2. 

Traverse the array from two directions using two index variables a and /3, 
initialized to 3 and n. The variable a increases until a point in P2 is observed; 

then (3 starts decreasing until a point in Pi is observed. Now, P[a] and P[l3] 
are swapped. The same process continues until a < j3. 

We use two scalar locations ni and n2 to store |Pi| and |P2|. We now sort both 
the set of points Pi and P2 separately with respect to their distances from iij. 
Note that, the distance values are not stored. While comparing a pair of points, 
their distance values are computed for the comparison. We now describe the 
method of computing all the PMERs with the points in Pi, keeping {pi,Pj) at 
its top boundary. 

As in the earlier algorithm, we consider a curtain whose two sides are bounded 
by the boundary of TZ, and top boundary is attached to both pi and pj. The 
curtain is allowed to fall. As soon as it hits a point p e Pi it reports a PMER. 
This point can easily be obtained from the sorted list of Pi, stored in the array 
P. If the projection tt of the point p on i^j lies inside the interval \pi,pj], the 
processing of iij stops. Otherwise, the curtain is truncated at tt, and the process 
continues to process the next point in Pi. After finishing the processing of all 
the points in Pi, we process the points in P2 in a similar manner to generate the 
PMERs with their bottom boundary passing through pi and pj. 



5.2 Correctness and complexity analysis 

The correctness of the algorithm follows from the fact that for each pair of points 
Pi,Pj & P, we have considered all possible PMERs with {pi,Pj) on its one side, 
and wc have considered each pair of points in P. 

Regarding the time complexity, note that, we have considered O(n^) pairs of 
points. For each pair, we have executed the procedure Process. Each time after 
the execution of the procedure Process, wc need to sort the array P with respect 
to their x and y coordinates for choosing the next pair of points for processing. 

In the procedure Process, the splitting of P into Pi and P2 needs 0{n) time. 
Sorting the members of Pi and P2 with respect to their distances from £ij needs 

O(nlogn) time, and then the reporting of the PMERs need 0(n) time. 

Note that, wc have only used four index variables z, 7, a and /3, and two integer 
locations ni and 77.2 to store size of Pi and P2- Thus we have the following 
theorem stating the complexity results of our proposed algorithm. 

Theorem 2. Given a set ofn points, the maximum area/perimeter rectangle of 
arbitrary orientation can be computed in 0{n^\ogn) time with 0(1) extra space. 
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